Molecular Recognition Effects in the Surface Diffusion of Large Organic Molecules: 

The Case of Violet Lander 
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Violet Lander (VL) (C108H104) is a large organic molecule that when deposited on Cu (110) 
exhibited lock-and-key like behavior (Otero et al, Nature Mater. 3, 779 (2004)). In this work we 
report on a detailed fully atomistic molecular dynamics study of this phenomenon. Our results show 
that it has its physical basis in the interplay of the molecular hydrogens and the Cu(110) atomic 
spacing, which is a direct consequence of an accidental commensurability between molecule and 
surface dimensions. This knowledge could be used to engineer new molecules capable of displaying 
lock-and-key behavior with new potential applications in nanotechology. 

PACS numbers: 68.65.-k, 61.46.-fw, 68.37.Lp, 71.15.-m 
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With the advent of nanoscience and nanotechnology 
and the perspective of molecular electronics [l|, H, H, 0, 
Hi H , significant theoretical and experimental efforts 
have been devoted to the study of the complex interac- 
tions involving organic molecular structures and metal- 
lic surfaces |& !, El d E E H \M- 0ne essential 
aspect of these phenomena is to understand how these 
interactions alter the properties of both molecule and 
surface. Recent progress has been achieved through the 
use of UHV-STM (ultrahigh vacuum-scanning transmis- 
sion microscopy) HH that allowed the identification 
of important structural and dynamical features related 
to the behavior of molecular wires adsorved on metallic 
surfaces. 

One important family of molecular wires is the so- 
called "Lander molecules" (because of its resemblance to 
a Mars surface rover) [llj • These large organic molecules 
are composed of a rigid polyaromatic ir central board and 
four spacers (legs) of up to eight 3,5-di-tert-butylphenyl 
cr-bonded to the central board (Fig. [l}. These spacers 
generate a configuration where the phenyl groups are 
nearly perpendicular to the main board plane by steric 
crowding. When deposited onto a metallic surface, these 
conformations allow an electronic decoupling (7r-bondcd 
from the metallic surface) to occur. Also, the presence of 
the tert-butyl groups, which increases the board-surface 
distance, permits some rotation of the "legs" without sig- 
nificantly reducing this distance. 

When adsorbed on Cu(100) or Cu(110) surfaces, 
the Lander molecule can act as a template for self- 
accomodating metal atoms at the step edges of the cop- 
per substrate, spontaneously generating metallic nanos- 
tructures dimensionally commensurable with the Lander 
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FIG. 1: Structural representation of the Violet Lander 
(C108H104) along different direction views. The experimen- 
tal and theoretical values of the major molecular dimensions 
L and W are displayed in Table I (see text for discussions). 



(la . Il8l | . This phenomenon opens up interesting possibil- 
ities for the bottom-up design of metallic nanostructures. 

Similar studies for other molecules of the Lander fam- 
ily, the so-called Violet Lander (VL) (Violet due to its 
color) - C108H104 - led to the observation in the solid state 
of the first non-biological lock-and-key like effect fl|| . A 
change by two orders of magnitude of the diffusion coef- 
ficients was observed when the molecular orientation of 
the substrate was changed. Previous experimental and 
theoretical investigations [HI, [2(| have revealed that 
VL molecules adsorb at room temperature (RT) on the 
Cu(llO) surface with the polyaromatic board parallel to 
the substrate and aligned along the [110] direction (par- 
allel configuration) [Fig. f^a), (c)]. In this configuration 
the molecules do not diffuse during the STM observation 
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FIG. 2: (color online) STM snapshots from the scanning of 
VL molecules adsorbed on Cu(llO) surface, (a)-(b) Molecules 
with their main boards aligned with the [110] direction (la- 
beled 1) do not diffuse; molecules rotated by 70° (labeled 
2) are able to diffuse along the [110] direction, (c)-(d) 3D- 
graphical atomistic representation of the VL in its (c) aligned 
and (d) 70° rotated configurations with respect to the [110] 
Cu(110) surface direction. 



time (several minutes). The estimated diffusion coeffi- 
cient is less than 5 x 10~ 19 cm 2 s _1 [l9j |. 

When the STM tip is used as a tool to push the ad- 
sorbed molecule onto the copper surface at low temper- 
atures_ (160-200 K) [lj], the molecule rotates 70° (from 
the [110] direction) [Fig. HJb), (d)], and the molecule dif- 
fuses along the [110] direction with an estimated diffusion 
coefficient of 4.8 x 10~ 17 cm 2 s _1 [lj|. This represents 
an increase of two orders of magnitude in relation to the 
parallel configuration. Interestingly, this molecular dif- 
fusion can be stopped at any time simply by flipping the 
molecule back to its parallel configuration. These two 
distinct configurations work as a two-state system (0 and 
1 / on and off), and it is the analogous of the so-called 
biological lock-and-key recognition between enzymes and 
the substrate on which they act [211 ]. 

Surprisingly, in spite of many years of investigation 
into the diffusion of organic molecules on metallic sur- 
faces, this phenomenon had not been observed before. 
As we shall see in the following discussions, its physical 
roots are based on an accidental commensurability of the 
VL and the atomic spacing of the Cu(110). Changing the 
molecule or/and the crystallographic direction, the phe- 
nomenon will not be observed. 

In this work we present a detailed fully atomistic 
molecular dynamics study of the diffusion process of VL 



TABLE I: Violet Lander dimensions (Fig. 1), in A, optimized 
with classical molecular mechanics (univeral force field 12211 ) . 
semi-empirical AMI method [pL DFT-LDA-Siesta [Mill, 
and DFT-LDA-Dmol3 H3,[^S. 



molecules adsorbed on a Cu(110) surface. We carried 
out molecular dynamics simulations in the framework of 
classical mechanics with standard force field [2^ |. which 
includes van der Waals, bond stretch, bond angle bend, 
and torsional rotation terms. This methodology has been 
proven to be very effective for the study of dynamical 
properties of complex structures [H, [24], H^, [26[ . 

For all simulations the following convergence criteria 
were applied: maximum force of 0.005 kcal/mol/A, root 
mean square (RMS) deviations of 0.001 kcal/mol/A, en- 
ergy differences of 0.0001 kcal/mol, maximum atomic dis- 
placement of 0.000 05 A, and RMS displacement of 0.000 
01 A. A selective microcanonical (constant number of 
particles, volume, and total energy) impulse dynamics 
was used, with time steps of 1 fs. 

Initially, the VL molecule was optimized in the gas 
phase (isolated). The use of full quantum methods for 
the whole system is not possible due to its large size, 
but the calculations for the isolated molecules is possi- 
ble. In order to test the quality of the geometrical re- 
sults for the molecular structures we have carried out a 
series of calculations using different methods: the semi- 
empirical hamiltonian AMI available in the GAMESS 
package [27j|. and the ab initio density functional meth- 
ods Siestajll, [H[ and state-of-the-art DMOL3 [H [H| 
in the local density approximation. The obtained struc- 
tural dimensions are very consistent and in very good 
agreement with the estimated experimental data from 
STM experiments [HI, [l9|, HH, [Hj, indicating that the 
used molecular force field reproduces the VL geometrical 
features quite well. 

Then the molecule with its board parallel to the 
Cu(110) surface is placed (about ~ 4-6 A) above the sur- 
face and set free to interact with the surface. We have 
also considered the cases where a little vertical impulse 
was given to the molecules to move them towards the 
surface. This is an additional test in order to provide 
enough kinetic energy to probe to local minima. These 
procedures were repeated varying the relative angular 
orientation and the geometries reoptimized. As the cop- 
per surfaces do not reconstruct, the Cu atoms were kept 
frozen at the experimental lattice value of a — 3.61 A 

Our results showed that the molecules always converge 
to two possible configurations, and 70° with respect to 
the [110] orientation. In order to better understand the 
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problem, we mapped the energy configuration as a func- 
tion of the rotation angle. For comparative purposes this 
was carried out in two different ways: the molecule was 
placed at 0° and initially with the board frozen, then ro- 
tated in steps of 5° and the geometry optimized. The 
procedure was then repeated, keeping frozen only the 
central ring of the board. The results are displayed in 
Fig. [31 and they are basically similar for the two pro- 
cedures. We observed that in fact the two most stable 
configurations are at and 70° (in agreement with the 
experimental data 19]), with a difference in energy of 
0.75 kcal/mol. A third minimum was obtained for about 
30° (but not experimentally observed), but considering 
its depth and width it could be overcome thermally. 

We repeated the process (frozen central ring and steps 
of 5°) in the [100] and [111] Cu surface directions. The 
results are displayed in Fig. 4. As we can see from Fig. 4, 
in contrast to the [110] surface (Fig. 3), the energy peaks 
and vallyes are smaller and smoother. No well defined 
constrained configuration was obtained. We observed in 
the simulations that these valleys and peaks could be 
easily thermically overcome. As discused in the following 
this energy could easily be overcome thermally for the 
[100] and [111] surfaces, which prevents the lock-key like 
phenomenon observed for the [110] surface. 
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FIG. 3: Relative total energy profile of a VL molecule de- 
posited on a Cu(110) surface as a function of the angle be- 
tween the main axis of the molecule and the [110] direction 
(see Fig. 2). For each angle, the molecule is optimized with 
its board or central ring being frozen. 

The results show that the energy necessary for the 
Cu(110) surface diffusion in the non-rotated case is much 
higher than in the rotated case, and beyond the thermally 
available energy of the temperature at which the experi- 
ments are realized, while for the rotated case the energy 
necessary is thermically available. 

Similarly to the rotational analysis (Fig. [3]), we have 
also mapped the translational (keeping only the central 
board ring frozen) movement for the rotated and non- 
rotated configurations of the Cu[110] surface. Due to 
their energy profile mentioned above, the translational 
mappings do not provide new relevant information about 
the (100) and (111) surfaces 
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FIG. 4: Relative total energy profile of a VL molecule de- 
posited on Cu[100] and Cu[lll] surfaces as a function of the 
angle between the molecule main axis and the [010] and [111] 
direction, respectively. 



The molecules were moved in steps of 0.1 A along the 
[110] direction, and for each point the geometry is opti- 
mized. The results are displayed in Fig. [5l Again, we ob- 
served that the energies associated with the movements 
of non-rotated and rotated cases are quite different, being 
higher for the non-rotated case. From Fig. 5 it is possi- 
ble to explain why the molecules can easily diffuse in the 
rotated case. The reason for this being that we have a 
very low energy barrier for diffusion. However, although 
the barrier is much higher in the non-rotated case than 
in the rotated case, its absolute value is not high enough 
to prevent the overcoming of the potential energy barrier 
by kinetic energy. Thus, although the static analysis gen- 
erated helpful information, it did not provide a complete 
physical description of the diffusion processes as the ex- 
plicit inclusion of dynamical aspects is needed in order 
to explain the experimental data. Static procedure anal- 
ysis has been often used in the literature. Our present 
analysis show that this procedure has limitations that 
could compromise the derived conclusions in some cases. 
These aspects have not been properly addressed before in 
the literature and might be of importance for the surface 
science of large organic adsorbates. 

A deeper understanding of the molecular diffusion can 
be obtained from the analysis of the temporal evolution 
through molecular dynamics simulations. Information on 
the force profile and molecular conformational changes in 
femto second scale can be easily addressed (not possible 
in the experimental case). 

We simulated the molecular diffusion of the VL 
molecules deposited onto a Cu(110) surface by impulse 
dynamics, i. e., we attributed an initial velocity to the 
molecule and followed its time evolution. This initial 
impulse is applied to mimic thermal effects or/and tip 
"kick" . We have considered the cases of the non-rotated 
and rotated configurations. In each case different impulse 
velocities were used to determine the threshold values to 
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FIG. 5: Energy profile for the VL molecule displacement onto 
the Cu(llO) surface, along the [110] direction in the rotated 
and non-rotated geometries. 
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FIG. 6: Force profile, as a function of time, as a result of the 
interaction between the VL molecule in the (a) non-rotated 
and (b) rotated geometries, and the Cu(110) surface. Times 
corresponding to 1.8, 4.8, and 15 ps as a result of the initial 
impulse are shown by arrows. 



induce molecular movements along the [110] direction. 
For the non-rotated case we concluded that initial val- 
ues larger than 0.9A/ps (equivalent to a kinetic energy 
of 1.36 kcal/mol) are necessary to induce molecular diffu- 
sion. On the other hand, for the rotated case only values 
of 0.4A/ps (equivalent to 0.27 kcal/mol) are required. 

In Fig. [5] we present the force profiles as a function 
of time of the forces experienced by the molecule for the 
different configurations. The displayed data are for a sit- 
uation where the initial impulse velocities were set up 
to 0.8A/ps along the [110] direction (see complementary 
material videoOl [36j]). For the non-rotated configuration 
(Fig. [(Jl), in the first ~ 1.8 ps the force exerted on the 
molecule attains high values near to 20 pN in the opposite 
direction to the initial impulse. The force values decrease 
quickly (~ 5.0 pN at 5.0 ps and < 2.5 pN at 15 ps) as 
the initial translational energy is converted into vibra- 
tional and torsional molecular movements. The molecule 
oscillates back and forth around its initial position, but 
no diffusion (net displacement) is observed (see videoOl). 
The situation is quite different for the rotated case (Fig. 
|BJd). The initial forces are quickly attenuated, and the 
molecule diffuses easily (see videoOl). 

From the videoOl we can clearly see that for the non- 
rotated case the molecule oscillates back and forth with- 
out diffusing. This oscillatory behavior can also be seen 
in the root mean displacement (RDM) data (Fig. 7). 
The associated RDM diffusion coefficients [13] are 9.1 x 
10 -6 cm 2 s -1 and 5.6 x 10~ 4 cm 2 s -1 for the non-rotated 
and rotated cases, respectively. These diffusion coef- 
ficients are obtained from the analysis of the impulse 
molecular dynamics trajectory simulations where the ini- 
tial impulse kinetic energy is quickly redistributed to the 
torsion/vibration/deformation- lengths molecular modes. 
In the experiments after the initial "tip kick" the non- 
rotated molecule undergoes a orientational transition to 
the rotated configuration, and due to the available ther- 
mal bath it has enough kinetic energy to diffuse. Al- 
though the absolute values of the theoretical and experi- 
mental diffusion coefficients cannot be directly compared, 
their differences can. They are in excellent agreement 
showing the same two orders of magnitude differences 
(9.1 x 10~ 6 cm 2 s _1 and 5.6 x 10~ 4 cm 2 s _1 versus 5.0 x 
10 _19 cm 2 s _1 and 4.8 x 10 _17 cm 2 s _1 for the non-rotated 
and rotated case, theoretical and experimental values, re- 
spectively, i. e., a factor of 100 difference. 

The detailed analysis of the molecular conformational 
changes as a function of time allowed us to get a clear 
picture of the diffusional processes. The relative position 
of the hydrogen atoms of the molecular "legs" in relation 
to the copper (110) direction plays a fundamental role in 
defining whether diffusion is possible or not. 

The VL adsorbed in the Cu(110) substrate exposes 
eight hydrogen atoms at the bottom of the legs. In 
the non-rotated configuration, these H atoms fit per- 
fectly ( "locked" ) into the fourfold hollow sites of Cu(110) 
(Fig. [Hk) ■ Due to this perfect fitting, any tentative trans- 
lational energy is more likely to be converted into confer- 




FIG. 7: Root mean displacement (RMD) of rotated (dot 
pointed curve) and non-rotated (fill curve) VL. The diffu- 
sion coefficient associated from curves are 9.1 x 10 _6 cm 2 /s 
and 5.6 x 10~ 4 cm 2 /s for non-rotated and rotated VL, respec- 
tively. 
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FIG. 8: (color online) Schematic view of a VL molecule in 
its (a) non-rotated, and (b) 70° rotated configurations with 
respect to the [110] Cu(110) surface direction. 

mational changes than into net displacement (videoOl), 
thus blocking diffusion. For the rotated configuration, as 
the fitting is not as good (Fig. [5b) as in Fig. [5k., transla- 
tional energy is easily converted into kinetic energy, and 
diffusion is possible. 

It is important to stress that the configuration, in 



FIG. 9: (color online) Schematic view of the VL molecule 
for the rotated configuration of (a) [110]; (b) [100]; and (c) 
[111] directions. The different hydrogen atom orientations in 
relation to the Cu surface atoms are clearly visible. 



which diffusion is blocked, is just a consequence of the 
accidental commensurability between the distance of the 
"legs" and the atomic spacing of Cu(110). For instance, 
for the Cu(lll) substrate the spacing is no longer com- 
mensurate with the "legs" and the perfect "fitting" no 
longer occurs, and consequently the lock-and-key like ef- 
fect ceases to exist (Fig. 9). The existence of the "lucky" 
commensurability in the case of VL/Cu(110) might ex- 
plain why this effect was not observed erlier. 

In principle, wc can engineer other molecules to satisfy 
the commensurability criteria to obtain other lock-and- 
key systems. As these molecules prove to be capable of 
exhibiting many interesting properties (such as sponta- 
neously building ordered nanostructures) , this observa- 
tion could be of great relevance to build nanostructures 
in a bottom-up approach. Wc hope the present study 
will stimulate further studies along these lines. 
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